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Abstract 

Two simple (rotator and one-particle) mechanistic models are suggested 
to describe simultaneously at a minimal level of sophistication two basic func- 
tions of Fi-ATPase: a motor regime driven by ATP hydrolysis and its inverted 
function as ATP synthesis. This description is consistent with the so-called 
rotary binding-change mechanism, a milestone of functioning ATP synthase, 
and uses a stepping (driving) function associated with two sequences of time 
instants, at which hydrolysis and synthesis reactions occur. It is useful to 
analyse experimental data and numerical simulations indeed predict corre- 
sponding dynamic behavior. 
PACS numbers: 05.60.-k, 05.40.-a, 87.10.+e 

Recently, the modeling of molecular motors - enzymes which transduce chemical energy 
into directed mechanical motion - based on the idea of gaining useful work by rectification 
of zero-mean noise fluctuations has attracted considerable attention JTJ and several models 
exploiting the so-called ratchet mechanism have been elaborated (see, e.g., Refs. PH^], to 



mention a few). In this context, ATP (adenosinetriphosphate) synthase shown schematically 
in Fiq. 1, being a realistic molecular engine of great interest nowadays [|Kj, should also be 
studied from the point of view of biological physics |TT]. This machinery is composed of 
two rotary motors: a membrane-embedded unit Fo and water-surrounded Fi-ATPase (also 
called Fi) connected by a coiled-coil 7 subunit. ATP synthase works as a reversible motor- 
pump machine: the proton flow through Fo is believed to create a clockwise (when Fi is 
viewed from the F side) torque that turns 7, whereupon ATP molecules are sequentially 
synthesized at three catalytic sites, one on each (3 subunit. Vice versa, ATP hydrolysis, 
which also occurs sequentially on /3's, but in the opposite direction, has been demonstrated 
to make 7 rotate backwards, converting F into a proton pump. In this case, 7 is driven 
by sequential hinge-bending motions of /3's, like a crankshaft by the pistons in a car engine 
(for more details see Ref. []12| ). In this paper, we focus on Fx-ATPase and consider its 
operation in both hydrolysis (as a motor) and synthesis (as a synthesizer) directions, using 
in parallel two simple mechanistic analogs: a plane rotator and a particle subjected to a 
periodic potential. Our description is consistent with Boyer's binding-change mechanism 
JTB|Ji"l|] , the recent findings of Cherepanov et al. and Junge et al. on elasticity properties 



of the 7 subunit, as well as with the experimental results of Yasuda et al. 

The structure of Fi-ATPase and the mechanics of motions within it are too complex to 
allow a detailed description of all interactions and motions of its different parts, i.e., the 



three /3 and one 7 subunits [p!0tp!2}^Ji7| --[19| . Therefore, it would be useful to describe this 
very sophisticated three-dimensional system, using simple mechanical equivalents (springs, 
particles, etc.) and keeping the main features of rotary dynamics found in previous studies, 



such as the modeling of Oster and Wang [12] and others []T0|,p~5|^19f| . This approach, "from 



complexity to simplicity", is often used in biological physics. The typical example of such 
a modeling is the propagation of a nerve impulse on the giant axon of the squid, where 
insight in the dynamics of the original Ho dgkin- Huxley equations [2(J has been obtained by 
reduction to the minimal FitzHugh-Nagumo system |?T] . 



Here the three-dimensional and four-body interaction in the ^37 subcomplex is suggested 



to be effectively described as a coupling of a planar rotator of length i? centered in the middle 
of an equilateral triangle (see left panels of Figs. 2 and 3) with three equivalent catalytic 
sites (denote them by numbers 1, 2, and 3) at the vertices of the triangle. Obviously, this 
rotor-stator interaction, resulting in a crankshaft-like rotation of 7 (and cooperative catalysis 
on /3's), is a periodic function: U{6 + 2n) = U{6), where 9 denotes the angular position of 
the rotator (positive if counter-clockwise). The mechanical equivalent of a driving torque 
on 7 can be chosen as a stretched spring or an effective particle displaced from equilibrium 
in the periodic potential U(6) as illustrated in the left and the right panels of Figs. 2 and 
3, respectively. At each site i — 1, 2, 3, a spring K{ is attached connecting this site with the 
rotator. All the three springs are supposed to be identical, but only one spring is allowed to 
be switched on, while the other two are switched off, at a time. Then sequential switching 
the springs will result in a power stroke on 7. Without loss of generality, the rotor-stator 
potential (given in units of the hydrolysis energy W ~ 80 pN-nm) can be written as 

U{6) = [r{6)-a?/{l-a)\ (1) 

where a = d/R is the length of each spring Ki being undistorted, r{6) = 



1 + (1 + a) 2 — 2(1 + a) cos 6 the instantaneous spring length, and I = r(±27r/3). 
Let us now describe how our spring system operates in both the hydrolysis and synthesis 
directions and how it can be coordinated with ATP hydrolysis/synthesis reactions at the 
catalytic sites of the j3 subunits. According to Boyer's hypothesis |L3j supported by the 



structural studies by Walker and coworkers |14|], and recently by direct observation of Noji 



et al. [|18j, each site 1, 2, or 3 can be found at least in one of the three states: T (ATP 



binding), E (empty), and D (ADP binding), at a time. Structurally, they are arranged as 
T, E, and D counter-clockwise (see Figs. 1-3) and can be put on the 9 axis at the lattice 
("catalytic") sites with spacing 2n/3, as shown in the right panels of Figs. 2 and 3. The 
dimensionless periodic (with period 2tt) potential U(6) [not necessary of the form ([I])] is 
supposed to be "rigidly tied" to this lattice, so that its minima are always found at sites 
with state T. We assume this potential to satisfy the normalization conditions: U(T)= 



and {7(D) = 1. The last constarint means that the potential energy of the effective particle 
in state D is equal to the free hydrolysis energy W. 

Consider first the hydrolysis direction when Fi operates as a motor (see Fig. 2). Let 
initially the 7 subunit (rotor) be found at equilibrium, performing there thermal fluctuations. 
This state is represented by spring K\ being undistorted and switched on state T, while 
the springs at the other sites {K2 and K3) are switched off. In the particle equivalent, this 
situation is represented in the right panel of Fig. 2(a) by a particle fluctuating in the vicinity 
of one of the minima of the potential {7(6?) , i.e., in state T. When an ATP molecule settles into 
site 2 which is found at this time in state E, it appears bound there, resulting in ATP binding: 



E — > T. According to the rotary binding-change mechanism |J13|]14| , this transition implies 
the two conformational changes at the next two sites: release of the inorganic phosphate 
(P) and the ADP (adenosinediphosphate) molecule from site 2 (D — > E), and hydrolysis of 
the ATP molecule located at site 3 (T — > D). After these state transitions have occured, 
spring K\ is switched off, while spring K 2 is switched on, causing the power stroke on the 
rotor as demonstrated by Fig. 2(b). Therefore the rotor is driven forward before it reaches 
a new equilibrium state, stepping forward by 27r/3. Correspondingly, as shown in the right 
panel, the potential U{6) steps forward (to the right) by 27r/3, so that the particle appears 
to be lifted uphill at the level {7 = 1, thereafter sliding down and finally dwelling in the 
next potential minimum before the next sequence of conformational changes takes place. 
However, if occasionally two sequential conformational transitions occur very close in time 
resulting in two almost simultaneous potential steps forward, the particle appears on the 
positive slope of the potential U(9), sliding thereafter down backwards. When the time 
between two sequential transitions is still short, but long enough for the particle to make 
a descent close to equilibrium, the double potential step forward will result in the double 
sliding down on the negative slope. Indeed, both such occasional steps of 7 loaded by an 
actin filament were observed in experiments flIS , [IB| , I9| . We denote the sequence of time 



instants when hydrolysis events occur by {tf}^ and call it a "hydrolysis sequence". 

Consider now the case, when an external torque (e.g., from the side of Fq) is applied to 



the 7 subunit clockwise, as shown in the left panel of Fig. 3(a). According to Cherepanov et 



al. and Junge et al. [OJ, the external energy of torsion is stored as an elastic strain energy of 
7. When the rotator turns clockwise totally by 2n/3, the torsional energy is liberated for the 
synthesis of one ATP molecule at the site being in state D. Then spring K% is switched off, 
while spring K 3 is on, bringing the system to the zero-energy level (U = 0). Equivalently, the 
particle "drops down" to equilibrium state T [see Fig. 3(b)]. These mechanical equivalents 
are again in accordance with Boyer fl3| for ATP synthesis, when the conformational change 



D — > T implies the transitions at the other two sites: T — > E and E — > D. We denote the 
sequence of time instants when synthesis events occur by {tj} c *L l and call it a "synthesis 
sequence" . 

Summarizing, one can present the diagram shown in Fig. 4, which describes the elastic 
strain energy of the torsional rotor-stator interaction as a function of 9, combining both the 
hydrolysis and synthesis directions. Arrows show the rotational direction under hydrolysis 
(9 > 0, motor regime) and synthesis (8 < 0, synthesizer regime) . 

In experiments |TB|,|TB|,|TP|] , the 7 subunit was loaded by an actin filament rotating in a vis- 
cous solution. Therefore, we approach the overdamped limit and the equation of motion for 
the rotator driven by a periodic potential U{6) fluctuating stepwise forward and backwards 
by 27r/3 (due to hydrolysis power strokes and an external load torque %) reads 

Y8 = -Wd e U[9 - 2irS(t)/3] - %{9, t) + £(*), (2) 

where T is a viscous friction coefficient, U(9) is normalized by U(2im) = and U[2iT(n ± 
1/3)] = 1, n = 0, ±1, . . . , and is the Brownian torque with the auto-correlation function 
(£(£)£(£')) = 2TkBT5(t — £'), where ks denotes Boltzmann's constant and T is the absolute 
temperature. The stepping (driving) function 

tf(f) = £ e (*-tf)-Ee(*-*7)> (3) 

i=l j=l 

where @(t) = for t < and 0(t) = 1 for t > 0, is defined through the two time sequences 
{tf}^ and {t~}^ =1 indicating the time instants when hydrolysis and synthesis reactions 
occur, respectively. 



The hydrolysis sequence {tf} is a random process determined by ATP concentration 
([ATP]). We define it through two characteristic duration times At D and At T by the recur- 
rence relation 

t+ x = t+ + At D + 2( i+1 At T , 2 = 1,2,..., (4) 

where each Q G [0, 1] is a random value with uniform distribution (since ATP concentration 
is supposed to be constant in the solution). The interval At T is large for low [ATP], but 
tends to zero as [ATP] is sufficiently high. Therefore one can assume that At T = A T /[ATP] 
with some constant A T > 0. The interval At D triggers the hydrolysis reactions by release of 
P and ADP on the next site of the stator. In the limit of low concentration of nucleotides 
(ATP and ADP), this interval is short, whereas for high concentration, the release of the 
hydrolysis products is impeded and therefore some saturation for At D takes place. As a 
result, one can assume that At D = At s t[ATP]/(C D + [ATP]) with some constant C D > 
and At s t being the duration time of one step. At zero temperature, it follows from the over- 
damped dynamics governed by Eq. (Q) that At s t — > oo (an overdamped particle approaches 
a potential minimum after an infinitely long time). However, at nonzero temperature, the 
time for one step becomes finite, because in the vicinity of equilibrium, the particle is "cap- 
tured" by thermal fluctuations. Thus, in the limit a — > 0, the potential ([p) is reduced to 
the simple form U = 2(1 — cos#)/3, admitting an explicit solution of Eq. (||]) for each step 
if £(£) = 0. Since the average amplitude of thermal fluctuations is J3V/2, one finds that 
At st ~ (3t /2) In (v / 3coty / 3P/8), where t = T/W is the time unit and V = k B T/W the 
dimensionless strength of white noise in Eq. (0). Since, at room temperature k B T ~ 4 
pN-nm, we thus have T> ~ 0.05, so that At s t can be estimated for each viscous load T. 

Averaging Eq. (||), one finds that {t^ +1 — tj) = At D + At T and therefore the rotational 
rate of 7 (number of revolutions per second) is V — [3 (At D + At x )] -1 . Inserting here the 
dependences of At DjT on [ATP] given above, one finds 

Vmax[ATP] 

K M +[ATP] 2 /(C D +[ATP]) ' 1 ) 

6 



where V max = (SA^)- 1 = (At D ) is a maximal average velocity (as [ATP] — > oo) and 
the constant K M = A T /At s t can be identified as the Michaelis constant, because Eq. ([5]) 
is reduced to the Michaelis-Menten law in the limit C D — > 0. According to Yasuda et al. 
fi~6|| , K M = 0.8 fiM and V max = 4 s -1 , and for these values the dependence (§) is shown in 
Fig. 5 (including the experimental data), with monotonic behavior for C D < K M . Similarly 
TBI , one can assume that the one-step duration depends on the length of the filament L as 



At s t — go + g±L 3 , and substituting this expression into Eq. (|), one finds for small C D the 
dependence V(L) ~ | (A T /[ATP] + g + giL 3 ) 1 plotted in Fig. 6, where the constants A T 
and g ,i are fitted to the experiments [16| . 

As described above, for a given random hydrolysis sequence (§), each synthesis instant 
tj is defined as the time when the rotor, being at some time in state T, rotates backwards 
by 27r/3. Direct simulations of the dimensionless (r = t/to) equation (0) with the potential 
([[]) are presented in Fig. 7. Here the intervals Ar D T = At D T /t in the sequence (f|) are given 
through the constants C D and K M , as well as the parameter Ar s t = At s t/t = W/3TV maiX to 
be evaluated from experiments. Thus, using that T(L = 1 /im) ~ 1 pN-nm-s and Knax = 4 
s _1 P^| , one finds Ar s t — 6.75. In the case without load (curves 1 and 2), the average 
velocity (9)/2itt = TV/W is in good agreement with the direct observations (see Fig. 3 of 
Ref. |16|) and the law plotted in Fig. 5 at C D = K M , where the two velocities shown with 
dotted lines 1 and 2 correspond to curves 1 and 2 in Fig. 7, respectively. Note that the 
constant C D in Eq. @ implicitly describes the main feature of trajectories at low [ATP] (see 
curve 2): on average the number of multi-steps exceeds that of steps backwards, as observed 
experimentally |18 , |16| . Next, when a load torque % > is applied, the rotational rate 
decreases with increase of this torque, as illustrated by curve 3 in Fig. 7. Moreover, when 
the load exceeds some threshold value, the motor operates in inverse, as shown by curves 4 
and 5. Similarly to Lattanzi and Maritan j8|, the law (|5|) can be modified by subtracting a 
positive constant that controls the direction of rotation if % > 0. 

Thus, we have developed two (rotator and one-particle) physical models of archetypal 



7 



simplicity, which are consistent with the rotary binding-change mechanism [|13],[14] and the 
elasticity properties of the 7 subunit fO^] . The cooperative rotational catalysis at the three 
j3 subunit s is described through two time sequences, each for switching hydrolysis and syn- 
thesis reactions, by adjusting the statistics of switching to satisfy the recent experimental 
results [frS|,rS . The models described in this paper are generic and simple; they do not 
depend on details of the periodic potential U(8). In the hydrolysis (motor) direction, the 7 
subunit works in a "passive" regime; only all the (3 subunits are coordinated in the cooper- 
ative rotational catalysis. In the reversible (synthesis) direction, 7 is "active", causing the 
corresponding (again cooperative, but in the inverse sequence) conformational changes after 
its strain energy of torsion reaches the free ATP hydrolysis energy. These important fea- 



tures are consistent with both Boyer's binding-change mechanism [0,|1J| and the findings of 
Cherepanov et al. and Junge et al. [|15| . The dependences of the model parameters on ATP 
concentration are general and physically motivated. In the framework of our description, 
the load torque Ti(9,t) generated by the F part can further be involved explicitly resulting 
in a general motor/pump model of ATP synthase. 

We also conclude that the puzzle, how does the binding-change mechanism work, may 
be essential not only for understanding the chemistry (dissipative catalysis) of creation of 
ATP, one of the most important processes in life, but also constitute a key physical problem 
behind the function of molecular motors, such as design of man-made molecular devices. 
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FIGURE CAPTIONS 



FIG. 1. Schematics of ATP synthase adapted from Ref. [|T^]. An asymmetric 7 shaft 
rotates relatively to the hexamer formed by a and three (3 subunits arranged alternatively. 
The other subunits, which constitute Fo including the "anchor" part are not shown. The 
positive direction of 7 rotation and the directions of proton flow and rotational catalysis 
(sequential synthesis/hydrolysis reactions in the hexamer) are shown by the arrows. 

FIG. 2. Spring (left panels) and particle (right panels) equivalents of the rotor-stator 
system evolving in hydrolysis direction, (a) Rotor (left) and particle (right) are found 
in equilibrium (in state T). (b) Power stroke caused by stretched spring K-i will rotate 7 
counter-clockwise (left). Sliding the particle down on the negative slope of periodic potential 
U(6), after it has moved forward by 27r/3 (right). 

FIG. 3. Spring (left panels) and particle (right panels) equivalents of the rotor-stator 
system evolving in synthesis direction, (a) External torque drives the rotator clockwise and 
elastic strain energy of the system is stored in spring K\ (left). Lifting particle uphill in 
potential U(9) (right), (b) Release of elastic strain energy after synthesis takes place (left). 
After the particle energy has reached the value {7=1, potential U(9) steps backwards by 
27r/3, allowing the particle to "drop downhill" to zero energy level (right). 

FIG. 4. Strain energy of the rotor-stator system against angular position 9 in hydrolysis 
(9 > 0) and synthesis (9 < 0). 

FIG. 5. Rotational rate of 7 against ATP concentration calculated for three values of 
C D , using Eq. fl5|). Experimental results [16| shown by circles are given for comparison. 

FIG. 6. Rotational rate of 7 against the length of actin filament: experimental results 
(circles, squares, and triangles) and dependences V(L) calculated for three values of ATP 
concentration. 

FIG. 7. Typical trajectories for different regimes of Fi. Curves 1 and 2 demonstrate pure 
motor regime (without load) for high (curve 1; At d = 6.75 and Ar T = 0) and low (curve 2; 
At d = 2, Ar T = 10, [ATP] = 0.4 fM, and C D = K M = 0.8 /iM) ATP concentration. The 
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other trajectories illustrate mixed motor /synthesizer regimes for different constant loads: 
below threshold (curve 3; At d = 1.7, At t = 20, and % = 20 pN-nm), nearby threshold 
(curve 4; At d = 6.75, Ar T = 0, and % = 7^ ~ 38.4 pN-nm), and above threshold (curve 5; 
At d = 6.75, Ar T = 0, and = 0.48 pN-nm). 
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FIG. 1. "Rotary Motion in F r ATPase: 
by A.V. Zolotaryuk et al. 




Fig. 2. "Rotary motion in ^-ATPase:...", by A. V. Zolotaryukef a/. 




Fig. 3. "Rotary motion in ^-ATPase:...", by A. V. Zolotaryukef a/. 




Fig. 4. "Rotary motion in ^-ATPase:...", by A. V. Zolotaryukef a/. 
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FIG.5. "Rotary Motion in F r ATPase: 
by A.V. Zolotaryuk et al. 
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FIG.6. "Rotary Motion in F r ATPase: 
by A.V. Zolotaryuk et al. 




Fig. 7. "Rotary motion in ^-ATPase:...", by A. V. Zolotaryukef a/. 



